AFRL-RZ-WP-TP-2010-2006 


STATOR  DIFFUSION  ENHANCEMENT  USING  A  RE¬ 
CIRCULATING  CO-FLOWING  STEADY  JET 
(POSTPRINT) 

David  Car,  Nicholas  J.  Kuprowicz,  Jordi  Estevadeordal,  Gecheng  Zha, 
and  William  Copenhaver 

Engine  Integration  and  Assessment  Branch 
Turbine  Engine  Division 


JUNE  2004 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 


STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
PROPULSION  DIRECTORATE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7251 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


NOTICE  AND  SIGNATURE  PAGE 


Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for  any 
purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U.S.  Government. 
The  fact  that  the  Government  formulated  or  supplied  the  drawings,  specifications,  or  other  data 
does  not  license  the  holder  or  any  other  person  or  corporation;  or  convey  any  rights  or 
pennission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  relate  to  them. 

This  report  was  cleared  for  public  release  by  the  Air  Force  Research  Laboratory  Public  Affairs 
(AFRL/PA)  Office  and  is  available  to  the  general  public,  including  foreign  nationals.  Copies 
may  be  obtained  from  the  Defense  Technical  Infonnation  Center  (DTIC)  (http://www.dtic.mil). 

AFRL-RZ-WP-TP-20 10-2006  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR  PUBLICATION 
IN  ACCORDANCE  WITH  THE  ASSIGNED  DISTRIBUTION  STATEMENT. 


*//Signature// 

NICHOLAS  J.  KUPROWICZ 
In-House  Program  Manager 
Engine  Integration  &  Assessment  Branch 
Turbine  Engine  Division 
Propulsion  Directorate 


//Signature// 

TIMOTHY  J.  LEWIS 
Chief 

Engine  Integration  &  Assessment  Branch 
Turbine  Engine  Division 
Propulsion  Directorate 


//Signature// 

ROBERT  D.  HANCOCK 
Acting  Principal  Scientist 
Turbine  Engine  Division 
Propulsion  Directorate 


This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange,  and  its 
publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its  ideas  or  findings. 

*Disseminated  copies  will  show  “//Signature//”  stamped  or  typed  above  the  signature  blocks. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 

Conference  Paper  Postprint 


1.  REPORT  DATE  (DD-MM-YY) 

June  2004 


4.  TITLE  AND  SUBTITLE 

STATOR  DIFFUSION  ENHANCEMENT  USING  A  RE-CIRCULATING  CO- 
FLOWING  STEADY  JET  (POSTPRINT) 


6.  AUTHOR(S) 

David  Car  (AFRL/RZTF) 

Nicholas  J.  Kuprowicz  (AFRL/RZTA) 

Jordi  Estevadeordal  (Innovative  Scientific  Solutions,  Inc.) 
Gecheng  Zha  (University  of  Miami) 

William  Copenhaver  (AFRL/RZT) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Fan  Compressor  Branch  (AFRL/PRTF)  Innova 

Engine  Integration  and  Assessment  Branch  (AFRL/RZTA)  2766  Ii 
Turbine  Engine  Division  (AFRL/RZT)  Beaver 

Air  Force  Research  Laboratory,  Propulsion  Directorate  — — 

Wright-Patterson  Air  Force  Base,  OH  45433-7251  niver 

Air  Force  Materiel  Command,  United  States  Air  Force  ..  .  , 


3.  DATES  COVERED  (From  -  To) 

01  April  2003  -  01  April  2004 


5a.  CONTRACT  NUMBER 

In-house 

5b.  GRANT  NUMBER 


Innovative  Scientific  Solutions,  Inc. 
2766  Indian  Ripple  Rd. 
Beavercreek,  OH  45440 


5c.  PROGRAM  ELEMENT 
NUMBER 

62203F _ 

5d.  PROJECT  NUMBER 

3066 _ 

5e.  TASK  NUMBER 

11 


5f.  WORK  UNIT  NUMBER 

30661117 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-RZ-WP-TP-20 10-2006 


University  of  Miami 

Dept,  of  Mechanical  Engineering 

Coral  Gables,  FL  33124^ 


10.  SPONSORING/MONITORING 

AGENCY  ACRONYM(S) 

AFRL/RZTA 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-RZ-WP-TP-20 10-2006 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 
Propulsion  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7251 
Air  Force  Materiel  Command 
United  States  Air  Force 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 

Conference  paper  submitted  to  the  Proceedings  of  the  ASME  Turbo  Expo  2004  Power  for  Land,  Sea  and  Air  held  in  Vienna,  Austria, 
June  14  -  17,  2004.  PA  Case  Number:  AFRL/PA  2003-3026,  25  Nov  2003.  Paper  contains  color.  The  U.S.  Government  is  joint 
author  of  this  work  and  has  the  right  to  use,  modify,  reproduce,  release,  perform,  display,  or  disclose  the  work. 


14.  ABSTRACT 

This  paper  outlines  a  steady  flow  control  technique  that  augments  the  diffusion  process  within  a  stator  passage  via  a 
continuous  co-flowing  secondary  flow  stream  along  the  suction  surface.  The  technique  is  similar  to  that  used  for  flow 
vectoring  in  nozzles  where  a  secondary  flow  stream  is  used  to  enhance  the  diffusion  and  vectoring  of  high  speed  jets. 
Diffusion  factors  in  excess  of  0.95  are  simulated  and  the  “penalty”  for  the  secondary  system  is  addressed  with  an 
availability  and  simple  power  analysis.  Losses  within  the  secondary  flow  stream  were  included  in  the  availability 
analysis,  but  it  did  not  account  for  losses  within  a  delivery  system  of  this  secondary  flow.  This  was  accomplished 
through  the  ID  power  analysis  which  assessed  this  technique’s  impact  on  the  efficiency  of  an  axial  compression  stage 
and  the  sensitivity  of  this  efficiency  to  the  secondary  flow  system’s  efficiency.  A  system  level  analysis  is  presented  to 
assess  the  specific  fuel  consumption  and  thrust-to-weight  merits  that  may  be  realized  in  a  notional  engine  with  this  type 
of  flow  control.  A  cascade  experiment  was  also  performed  to  demonstrate  the  concept  in  a  blow-down  cascade  tunnel. 

15.  SUBJECT  TERMS 

flow  control,  stator  diffusion,  compression,  turbine  engine 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

17.  LIMITATION 
OF  ABSTRACT: 

SAR 


18.  NUMBER 
OF  PAGES 

20 


19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

Nicholas  J.  Kuprowicz 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 
N/A 


1 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


Proceedings  of  ASME: 
ASME  Turbo  Expo  2004  Power  for  Land,  Sea  and  Air 
June  14-17,  2004-  Vienna,  Austria 


GT2004-53086 


STATOR  DIFFUSION  ENHANCEMENT  USING  A  RE-CIRCULATING  CO-FLOWING  STEADY  JET 


David  Car 
AFRL/PRTF 
1950  Fifth  St.  Bldg.  18D 
Wright-Patterson  AFB,  OFI  45433 


Nicholas  J.  Kuprowicz 

AFRL/PRTA 
1950  Fifth  St.  Bldg.  18D 
Wright-Patterson  AFB,  OH  45433 


Jordi  Estevadeordal 

Innovative  Scientific  Solutions,  Inc. 
2766  Indian  Ripple  Rd. 
Beavercreek,  OH  45440 


Gecheng  Zha 

University  of  Miami 
Dept,  of  Mechanical  Engineering 
Coral  Gables,  FL  33124 


William  Copenhaver 
AFRL/PRT 

1950  Fifth  St.  Bldg.  18D 
Wright-Patterson  AFB,  OH  45433 


ABSTRACT 

This  paper  will  outline  a  steady  flow  control  technique  that 
augments  the  diffusion  process  within  a  stator  passage  via  a 
continuous  co-flowing  secondary  flow  stream  along  the  suction 
surface.  The  technique  is  similar  to  that  used  for  flow 
vectoring  in  nozzles  where  a  secondary  flow  stream  is  used  to 
enhance  the  diffusion  and  vectoring  of  high  speed  jets. 
Diffusion  factors  in  excess  of  0.95  are  simulated  and  the 
“penalty”  for  the  secondary  system  is  addressed  with  an 
availability  and  simple  power  analysis.  Losses  within  the 
secondary  flow  stream  were  included  in  the  availability 
analysis,  but  it  did  not  account  for  losses  within  a  delivery 
system  of  this  secondary  flow.  This  was  accomplished  through 
the  ID  power  analysis  which  assessed  this  technique’s  impact 
on  the  efficiency  of  an  axial  compression  stage  and  the 
sensitivity  of  this  efficiency  to  the  secondary  flow  system’s 
efficiency.  Also,  a  system  level  analysis  is  presented  to  assess 
the  merits  that  may  be  realized  in  a  notional  engine  with  this 
type  of  flow  control.  Particularly,  impacts  on  specific  fuel 
consumption  and  thrust-to-weight  ratio  were  addressed.  A 
cascade  experiment  was  performed  to  demonstrate  the  concept 
and  was  conducted  in  a  blow-down  cascade  tunnel.  Significant 
improvements  in  diffusion  were  qualitatively  seen  from  the 
DPIV  measurements  despite  limitations  in  achieving  the 
desired  secondary  flow  conditions. 

INTRODUCTION 

The  trend  to  higher  stage  loading  in  axial  turbomachines 
continues  to  put  an  increasing  demand  on  the  stator  to  provide 
higher  amounts  of  diffusion.  Many  researchers  are  addressing 
this  issue  through  the  use  of  varying  flow  control  techniques  [1- 
4],  Many  of  these  techniques  involve  the  use  of  blowing  jets  or 


suction  devices  or  a  combination  of  the  two  to  control  the 
amount  of  separation  and  thereby  increase  the  diffusion  and 
overall  performance  of  the  blade  section.  There  are  others  who 
are  showing  promise  using  unsteady  techniques  with  synthetic 
jets  to  control  boundary  layer  parameters  [5,6],  These 
techniques  are  very  appealing  in  that  they  impart  a  zero  net 
mass  change  to  the  core  flow.  Still  others  are  using  plasmas  for 
boundary  layer  control  [7]. 

The  flow  control  technique  outlined  in  this  paper  imparts 
no  net  mass  flow  change  to  the  core  flow  through  the  use  of  a 
re-circulating  co-flowing  jet  that  enhances  the  diffusion  and 
turning  levels  through  its  interaction  with  the  core  flow.  This 
technique  is  similar  to  those  used  in  jet  nozzle  flow  vectoring 
which  create  a  secondary  flow  using  a  suction  collar  and  pump 
[8-11],  The  technique  outlined  here  spawned  from  the  work 
performed  in  fluidic  vectoring  in  nozzles  and  is  an  attempt  to 
create  a  similar  condition  in  the  stator  passage  with  a  re¬ 
circulating  co-flowing  system.  It  is  envisioned  that  this 
secondary  flow  is  sustained  with  a  re-circulating  secondary 
flow  system  as  shown  in  Figure  2. 

The  “penalty”  of  this  type  of  arrangement  on  the  overall 
compression  system  may  not  be  too  great  if  the  required  total 
pressure  ratio  of  the  secondary  flow  system  and  its  efficiency 
requirement  is  reasonable.  A  power  and  availability  analysis 
was  conducted  in  an  attempt  to  address  this  issue.  Also,  a 
system  level  analysis  using  a  notional  engine  configuration  is 
carried  out  to  assess  the  impact  on  engine  performance.  The 
concern  over  heating  effects  caused  by  the  re-circulating  flow 
and  the  possible  need  for  intercooling  the  secondary  flow  was 
raised  in  private  communications  [12],  This  technique  relies  on 
momentum  and  heat  exchange  between  the  core  flow  and  the 
secondary  flow  as  it  traverses  the  blade  so  it  may  not  be 
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necessary  for  an  intercooler,  particularly  if  the  secondary  flow 
system  does  not  require  a  substantial  power  input. 

A  cascade  experiment  was  conducted  to  evaluate  the 
concept.  Shortfalls  within  the  secondary  flow  system  limited 
the  amount  of  secondary  flow,  falling  short  of  the  desired 
conditions  of  the  CFD  simulations,  but  a  qualitative 
improvement  in  diffusion  was  apparent  from  DPIV 
measurements. 


C,2 


CONTROL  VOLUME  ANALYSIS 

A  2D  control  volume  analysis  was  performed  to  determine 
the  relevant  non-dimensional  parameters  and  the  influence  this 
flow  control  technique  may  have  on  the  core  flow.  The  control 
volume  used  is  shown  in  Figure  3.  The  simplifying  assumption 


is  made  that  the  region  containing  the  secondary  flow  jet  is  a 
circular  arc  and  that  the  inlet  and  exit  velocity  of  the  secondary 
flow  jet  is  equal.  The  resulting  non-dimensional  axial  and 
tangential  blade  loading  equations  that  result  are  as  follows: 

K  =  1  -  -  CM  (sin(6»2 )  -  sin(^  ))  -  C 

Vxl  (1) 

F*y  =  r'  +  CM  (cos(#2 )-  cos(^, )) 

where 

c  _  PjVfh  _  riijVj 
M  PKS  mcVxl 

r.  (yyi-yy2)s_  r 
Vxls  Vxls 


c 


p 


P2~P1 

pK 


Figure  3  Control  Volume 

Albeit,  the  reader  will  recognize  these  parameters  as  the 
static  pressure  rise  coefficient,  C  ;  the  circulation,  T  ;  the 

momentum  coefficient;  C  ;  and  the  non-dimensional  tangential 

and  axial  blade  loading,  r  and  r  x  respectively. 

The  amount  of  turning  a  stator  cross  section  achieves  is 
proportional  to  the  tangential  blade  loading.  The  momentum 

coefficient,  C ^ ,  offers  a  positive  contribution  to  the  tangential 

* 

blade  loading,  r  ,  if  and  only  if  it  effects  1  in  a  positive 
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sr  n 

manner,  i.e.  - >  (J  .  Interestingly,  if  the  secondary  flow  jet 

dcp 

direction  were  reversed  in  the  control  volume  analysis,  i.e.  a 

counter-flow  direction,  then  the  equations  remain  the  same. 
_  * 

Also,  r  is  increased  if  the  angular  extent  over  which  the 
secondary  flow  region  acts  is  increased. 


The  momentum  coefficient,  c  ,  offers  a  scaleable  non- 

dimensional  parameter  for  the  secondary  flow  system.  It’s  seen 
in  much  literature  concerning  flow  control  devices;  particularly 
ones  involving  circulation  control  airfoils  [13,14],  The 
secondary  mass  flow  must  scale  with  the  core  mass  flow  along 
with  the  jet  velocity  to  axial  velocity  ratio  in  order  to  hold  the 
momentum  coefficient  constant.  The  secondary  mass  flow 
could  also  scale  with  the  engine  shaft  speed  since  the  core  mass 
flow  in  an  axial  compression  system  is  tied  to  the  engine  shaft 
speed  and  it’s  envisioned  that  the  secondary  flow  pump  would 
be  tied  to  the  same  shaft  as  shown  in  Figure  2.  Therefore  it 
seems  reasonable  that  the  type  of  configuration  shown  in  Figure 


2  could  hold  c  relative  constant  over  the  engine  operating 

range.  A  computational  fluid  dynamics  investigation  was 
performed  as  a  first  step  in  understanding  the  influence  this 
flow  control  technique  might  have  on  the  core  flow. 


CFD  SIMULATIONS 

The  geometry  shown  in  Figure  3  was  simulated  using 
FLUENT,  a  commercial  CFD  package.  The  geometry  was 
created  from  a  series  of  three  circular  arcs  ensuring  continuity 
of  slope  between  them.  The  first  circular  arc  extends  from  0- 
10%  axial  chord;  the  second  from  10-90%;  and  the  third  from 
90-100%.  The  overall  chord  is  0.127  m  with  a  solidity  of  1.67. 
The  radius  of  curvature  for  each  section  is  as  follows:  0.1836 
m,  0.0696  m  and  0.1841  m.  The  thickness  distribution  for  the 
controlled  blade  was  set  to  allow  enough  cross  sectional  area  to 
pass  the  required  secondary  flow  for  a  future  cascade 
experiment  which  would  allow  full  span  control.  This 
experiment  will  be  discussed  later  in  this  paper. 

The  simulation  was  a  2D  simulation  with  a  hybrid  grid 
consisting  of  triangular  elements  in  the  core  flow  with  a 
quadrilateral  element  O-grid  around  the  blade.  The  grid  system 
is  shown  in  Fig.  4.  The  solver  setup  used  k-epsilon  turbulence 
modeling  and  wall  functions  to  model  the  boundary  layer.  For 
the  flow  controlled  blade,  the  inset  cavity  grid  was  21x101  and 
there  were  a  total  of  17050  cells  in  the  entire  grid  domain. 

The  inlet  conditions  for  this  cascade  section  were  a  Mach 
number  of  0.7  with  an  inlet  angle  of  68  degrees.  An  axial 
exiting  flow  with  an  approximate  Mach  number  of  0.2  was  the 
desired  exit  condition.  This  case  would  result  in  an  extreme  DF 
of  approximately  0.98.  The  inlet  Mach  number  of  0.7  was 
chosen  because  of  a  future  cascade  test  in  which  the  cascade 
tunnel  was  limited  to  a  maximum  running  Mach  number  of  0.8. 


Figure  4  Hybrid  CFD  Grid 


The  inlet  and  exit  of  the  core  flow  domain  and  the  inset 
cavity  were  simulated  in  the  same  manner.  Both  specified  the 
total  temperature,  pressure  and  flow  angle  at  the  inlet  and  the 
static  pressure  at  the  exit  boundary.  The  exit  static  pressures  to 
the  core  flow  and  the  cavity  flow  were  iterated  until  the  desired 
inlet  Mach  number  to  the  core  flow  was  achieved  and  the  exit 
mass  flow  to  the  inset  cavity  was  equal  to  the  inlet  mass  flow  to 
the  cavity.  The  inlet  and  exit  to  the  inset  cavity  are  decoupled 
in  this  way  which  is  different  from  the  envisioned  application 
requiring  a  re-circulating  pumping  system.  Time  constraints 
forced  us  to  use  the  existing  boundary  conditions  provided  in 
the  CFD  package  rather  than  write  custom  routines  to  represent 
the  pumping  system  and  provide  feedback  to  the  inlet  of  the 
inset  cavity.  Therefore,  the  total  temperature  for  the  cavity 
inlet  was  set  to  the  core  flow  total  temperature  thereby  not 
addressing  the  possible  re-heating  effects  mentioned  earlier. 

First,  a  baseline  was  run  to  compare  against  and  establish 
the  effectiveness  of  the  flow  control  technique.  The  baseline 
blade  was  similar  in  shape  to  the  flow  control  blade  in  terms  of 
blade  curvature.  The  suction  surface  of  the  baseline  blade  did 
not  contain  the  inset  cavity  and  the  pressure  surface  was  altered 
on  the  flow  control  blade  to  provide  the  necessary  cross 
sectional  area  for  internal  cavities  for  the  suction  and  blowing 
volume  in  a  future  cascade  experiment. 

The  comparison  between  the  flow  controlled  blade  and  the 
baseline  blade  shape  is  still  valid  even  though  the  blade  shapes 
are  not  exactly  the  same  because:  1)  If  I  were  to  simulate  a  true 
baseline,  I  would  simulate  the  exact  flow  control  blade  with  the 
inset  cavity.  No  one  would  design  a  blade  shape  with  a 
backward  facing  step  on  the  suction  surface.  Therefore,  I 
would  have  to  simulate  a  blade  shape  without  the  inset  cavity 
which  changes  my  flow  controlled  blade;  2)  This  flow  control 
technique  inherently  modifies  the  blade  shape  by  needing  an 
inset  cavity,  unlike  other  techniques  that  seek  to  influence  the 
flow  by  applying  suction  and  blowing  to  an  already  existing 
profile.  The  technique  would  inextricably  alter  any  blade  shape 
it  was  applied  to.  With  reasons  1)  and  2)  stated  above,  the  goal 
should  be  to  compare  the  flow  control  technique  to  a  blade 
shape  that  tries  to  achieve  the  same  end.  Since  the  baseline 
shape  was  used  to  create  the  flow  controlled  shape,  particularly 
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keeping  the  curvature  reasonably  the  same,  the  comparison  is 
reasonable. 

The  baseline  case  was  simulated  and  the  Mach  contours 
are  shown  in  Figure  5.  As  can  be  expected,  the  flow  is  badly 
separated  on  the  suction  surface  starting  at  approximately  50% 
chord.  The  flow  controlled  blade  is  shown  in  Figure  6  for  a 
momentum  coefficient  of  0.197.  The  flow  remains  attached  for 
the  entire  blade  region  and  appears  completely  diffused  at  the 
exit  plane. 


Figure  5  Baseline  Blade  Mach  Contours 


Figure  6  Flow  Control  Blade  Mach  Contours 

Table  1  summarizes  the  baseline  and  flow  controlled  blade 
for  various  momentum  coefficients  (which  were  varied  by 
changing  the  secondary  flow  injection  to  core  flow  total 

pressure  ratios ,7Tjc )  and  fractional  injection  slot  heights  (h/S). 

Referencing  the  locations  highlighted  in  Figure  6,  the  definition 
of  DF  used  in  the  table  is: 


DF  =  1 


c,  3 


AV 


+  - 


V„ ,  2oV„ 


(2) 


where  ACV  is  the  change  in  tangential  velocity  between 

locations  (c,2)  and  (c,3).  Mass  averaged  quantities  were  used 
for  the  calculation. 


The  zero  momentum  coefficient  case  is  for  the  baseline 
solution.  As  the  table  shows,  the  circulation  and  hence  the 
blade  loading  is  increased  with  increasing  values  of 

injection  c  .  The  flow  controlled  blade  was  able  to  achieve  an 

axial  exit  condition  for  an  injection  c  between  0.179  and 

0.197  for  the  h/s=0. 011974  and  between  0.134  and  0.146  for 
h/s=0. 00661.  The  circulation  and  DF  values  in  Table  1  are 
plotted  in  Figure  7.  DF  will  be  discussed  as  a  figure  of  merit 
throughout  the  body  of  text  since  as  mentioned  in  the 
introduction,  the  desire  to  increase  axial  compression  stage 
loading  places  an  increased  demand  on  the  diffusion 
requirements  of  the  stator. 

Table  1  Flow  Control  Parameters 
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Figure  7  Circulation  and  DF  versus  Momentum  Coefficient 


The  slope  of  the  circulation  line  as  a  function  of  c  is 

dr*  A 

positive,  i.e.  - >  I) .  Therefore,  as  expected  the  influence 

dcM 

on  the  loading  is  beneficial  and  the  overall  diffusion  is  greatly 
increased.  The  relationship  is  also  seen  to  be  quite  linear  in  the 
range  simulated,  although  more  points  would  be  necessary  of 
course  to  establish  a  trend. 

The  DF  and  circulation  in  Figure  7  do  not  collapse  on  the 
injection  momentum  coefficient  for  the  two  different  injection 
slot  heights  simulated.  As  can  be  seen,  an  offset  exists  between 
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the  two  cases.  This  can  be  understood  by  realizing  that  the 
flow  control  process  does  not  rely  solely  on  injection,  but  also 
on  the  applied  suction  and  that  the  momentum  coefficients  for 
these  two  areas  are  quite  different.  The  preceding  control 
volume  analysis  assumed  that  the  inlet  and  exit  heights  (i.e. 
areas)  were  the  same  for  the  blowing  and  suction  ports  in  the 
inset  cavity.  In  reality  this  is  not  the  case.  The  suction  area 
needed  to  be  larger  than  the  inlet  area  to  prevent  choking  and 
accommodate  the  required  mass  flow  injected  and  mixed  with 
the  core  flow.  Therefore,  the  above  plot  is  particular  to  the 

geometry  at  hand  since  the  x-axis  is  c  for  the  injection  only 


and  does  not  take  into  account  the  separate  c  for  the  suction. 


Another  interesting  result  is  that  the  c  =0.146  with 

h/s=0.00661  achieves  a  similar  DF  as  the  0.197  with 
h/s=0. 01974  case.  Similar  results  were  discovered  in  the  past 
for  flow  control  airfoils  [13,14]  which  achieved  similar  lift 
coefficients  for  varying  injection  slot  heights  and  therefore 
various  momentum  coefficients.  It  was  a  matter  of  whether  you 
are  pressure  limited  or  flow  limited  in  the  secondary  flow 
system.  A  low  pressure  with  a  higher  flow  rate  could  achieve 
the  same  end  as  a  high  pressure,  low  flow  rate  system.  Of 
course  there  is  most  likely  a  limit  in  the  range  over  which  this 
applies  and  certainly  there  are  secondary  effects  that  would 
make  one  option  more  appealing  over  another,  but  this  appears 
to  be  case  for  the  narrow  range  of  geometries  and  momentum 
coefficients  simulated  here. 

A  parameter  that  helps  to  collapse  the  curves  is  the 
isentropic  power  input  of  the  secondary  flow  stream.  To  tie  in 
this  analysis  with  a  subsequent  stage  analysis,  Figure  2  and 
Figure  6  will  be  used  for  the  subscript  nomenclature.  The 
projection  of  the  locations  specified  in  Figure  2  to  the  CFD 
simulation  is  shown  in  Figure  6.  The  core  inlet  and  exit 
conditions  are  locations  (c,2)  and  (c,3)  respectively.  Likewise, 
(j,l)  and  (j,2)  represent  the  injection  and  suction  of  the 
secondary  flow  respectively.  The  fluidic  power  associated  with 
the  total  pressure  drop  between  the  suction  and  injection  ports 
can  be  defined  as: 


Pjjsen  =  ™jCpTtJ,  2 


^  y—\ 

7T  : 


-1 


(3) 


The  total  pressure  ratio,  71- ,  in  Eq.  (3)  is  the  ratio  of  the 
injection  to  suction  slot  total  pressure.  Non-dimensionalizing 
this  equation  using  the  inlet  enthalpy  flux,  rhcC  Tt  c  2 ,  yields: 

(  y- 1  \ 


P*  =/r. 

Kisen  J  i 


TC  ; 


-1 


(4) 


V  7 


where  T  -c  is  the  temperature  ratio  between  secondary  flow 

suction  slot  and  the  inlet  core  conditions  and  f  is  the  flow 
fraction,  i.e.  ratio  of  secondary  to  core  mass  flow. 

Plotting  the  DF  against  Eq.  (4)  using  mass  averaged  values 
is  shown  in  Figure  8.  As  is  seen,  the  DF  collapses  rather  well 
for  the  two  simulations  with  varying  injection  slot  height.  It 
should  be  noted  that  the  suction  slot  height  was  1.61  times 
larger  than  the  injection  slot  for  the  h/s=0. 00661  case  and  2.16 
times  larger  for  the  h/s=0.011974  case.  Interestingly,  this 


suggests  that  the  diffusion  enhancement  using  this  technique 
correlates  with  the  irreversibility  (power  loss)  of  the  secondary 
flow  stream.  Recently,  secondary  flow  stream  irreversibility 
has  been  presented  as  a  mechanism  for  the  flow  vectoring  using 
a  counter-flowing  secondary  flow  stream  for  nozzle  flow 
control  [15].  The  idea  is  presented  that  losses  within  the 
secondary  stream  due  to  the  interaction  with  the  core  flow 
increase  the  loading,  thereby  vectoring  the  flow. 

Looking  back  at  Table  1  again  also  shows  that  a  flow 

fraction  of  4.87%  is  required  for  a  C  of  0.197.  A  rule  of 

thumb  that  has  emerged  over  the  years  has  been  that  a  flow 
control  device  that  uses  0.5  %  of  the  core  flow  is  excessive. 
This  stems  from  the  understanding  that  0.5%  of  the  flow  would 
need  to  be  removed  or  added  elsewhere  in  the  cycle  and  this 
would  be  detrimental  on  the  overall  engine  performance.  But, 
unlike  other  techniques,  the  flow  is  re-circulated  rather  than 
added  or  removed.  In  order  to  obtain  a  better  understanding  of 
whether  this  amount  of  flow  fraction  would  cause  a  large 
penalty  on  the  compression  system,  both  an  availability 
analysis  and  a  simplified  power  analysis  were  performed. 
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Figure  8  Isentropic  Power  Input  versus  DF 


AVAILABILITY 

A  treatise  on  availability  analysis  is  given  in  [16].  The 
availability  of  a  system  is  defined  as: 

y¥  =  m(AHt-TtAs)  (5) 

where  A H  t  is  the  change  in  stagnation  enthalpy,  Tt  is 

initial  starting  total  temperature  and  As  is  the  change  in 
entropy.  The  availability  is  simply  the  available  power 
extractable  at  the  exit  to  a  system  if  the  system  were  brought 
isentropically  back  to  rest  at  its  initial  conditions.  For  a 

cascade  configuration,  generally  'P  <  0  since  the  cascade  is  a 
only  a  loss  producing  mechanism  with  no  power  input,  i.e. 
AHt  =  0  and  A.V  >  0  .  The  total  availability  of  a  system  is 
the  sum  of  the  availability  at  each  boundary  to  the  system. 
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As  before,  this  analysis  will  refer  to  Figure  2  and  Figure  6 
for  the  subscript  notation.  The  availability  for  the  CFD 
simulation  at  each  location  is  as  follows: 


w  = 

xc,  2 

0 

^,3  = 

»V 

= 

0 

VF. ,  = 

m . 

j 

Both 

c,2  £ 
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\Sc,3  -Sc.2 


)) 


(6) 


\SJ,2~Sj 


,)) 


these  two  streams  are  the  conditions  at  (c,2)  and  (j,l) 
respectively.  Summing  the  above  relationships  to  obtain  the 
total  availability,  simplifying  and  non-dimensionalizing  by  the 
core  inlet  enthalpy  flux  yields: 


'F 


rhccp,cTt,c,2 


T- 1  - 


As,, 


+  f 


cp.j  T>J, t 


CP,C  T,.c, 2 


As. 

r,  -1 - ^ 


-|(7) 


'p.f  7 


where  is  the  specific  heat  at  constant  pressure  and 


r 


c 


For  the  CFD  simulation,  the  secondary  flow  injection  total 
temperature  was  equal  to  the  core  flow  total  temperature,  i.e. 

Tt  j ,  =  Tt  c  2  .  Also,  for  simplification,  let  C p  c  =  C p  •  and  the 


Eq.  (7)  reduces  to: 

¥,=(rI-l)+/(r,-l)-^  +  /A,;-*  (8) 

Cp,c 

The  first  two  terms  in  Eq.  (8)  represent  the  work  input/loss 
of  the  core  and  secondary  stream  respectively  and  the  last  term 
the  combined  irreversibility.  The  work  occurs  through  the 
shear  work  between  the  core  and  the  secondary  flow  and 
although  small,  is  still  on  the  order  of  the  irreversibility  term 
[17],  Therefore,  it  must  be  included  in  the  analysis. 

The  ratio  of  availability  to  the  baseline  availability  (i.e.  the 
case  with  no  flow  control)  is  shown  in  Figure  9  for  two 
different  injection  slot  heights.  Mass  averaged  values  were 
used  for  the  calculations  and  the  exit  conditions  were  integrated 
at  the  constant  axial  location  labeled  (c,3)  in  Figure  6.  Note 
that  a  mass  averaged  value  of  the  baseline  case  in  this  region 
would  be  a  conservative  estimate  of  the  loss  since  the 
weighting  would  be  higher  for  the  higher  momentum,  i.e.  non- 
separated,  regions  of  the  flow.  In  other  words,  losses  would  be 
larger  is  this  flow  condition  were  allowed  to  mix  out. 
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Figure  9  Availability  versus  Momentum  Coefficient 


As  you  can  see  from  Figure  9,  for  the  highest  momentum 
coefficient  the  availability  of  the  flow  control  case  is 
approximately  50%  less  than  that  of  the  uncontrolled  case. 
This  simply  means  there  is  a  50%  reduction  in  the  overall  loss 
production  with  the  entrainment  and  mixing  process  involved 
for  this  flow  control  technique  as  compared  to  the  separated 
baseline  case.  Thermodynamically  this  flow  control  process  is 
more  efficient  than  the  separated  baseline  case,  but  losses 
within  a  delivery  system  could  make  the  two 
thermodynamically  equivalent.  The  next  section  will  evaluate 
the  performance  of  a  compression  stage  that  may  incorporate 
this  flow  control  scheme  and  assess  the  sensitivity  of  the  stage 
performance  to  the  secondary  flow  delivery  system  efficiency. 
Keep  in  mind  though,  even  it  were  determined  that  the  flow 
control  scheme  was  found  to  be  thermodynamically  equivalent 
with  the  separated  baseline  case  due  to  the  inclusion  of  losses 
within  the  delivery  system,  the  flow  controlled  case  still  has  the 
advantage  of  achieving  the  desired  diffusion  levels  and  exit 
conditions  over  the  uncontrolled  baseline. 

Another  consideration  from  Figure  9  is  that  the  power  loss 

for  the  lower  C  is  less  than  the  higher  C  values.  Since 


vjD 


'F; 


=  1  when 


base 


C  =  0  then  there  must  exist  a  value  of  C'u 


where  the  power  loss  is  a  minimum,  but  does  not  coincide  with 
the  value  of  C/(that  produces  the  most  diffusion  and  turning. 

This  may  be  an  artifact  of  the  averaging  procedure,  but  area 
averaging  produced  the  same  trend. 


POWER  ANALYSIS 

An  assessment  of  this  technique  used  in  the  stage  of  an 
axial  compression  system  was  performed  by  analyzing  the 
power  input  of  a  stage  using  this  flow  control  technique  in  the 
stator  section  as  shown  in  Figure  2  and  comparing  that  with  the 
isentropic  power  input  for  the  same  pressure  ratio.  To  obtain 
the  stage  design  characteristics,  a  meanline  analysis  was 
performed  with  the  goal  of  a  68  degree  rotor  exit  condition. 
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1.67  stator  solidity  and  0.98  stator  DF  in  order  to  be  consistent 
with  the  previous  results  of  the  CFD  analysis.  The  following 
meanline  input  parameters  were  used: 


Table  2  Meanline  Input  Parameters 


777 

195.3  kg/s/m- 

A 

Rotor  Tip  Speed 

426.72  m/s 

Rotor  Work  Input 

1.62820e5  J/kg 

AR  rotor 

1.1 

V  rotor 

0.90 

R  rotor 

0.85 

CJ  , 

rotor 

3.0 

CJ  ,  , 

stator 

1.67 

A  VR stator 

0.8 

The  resulting  pressure  ratio  and  reaction  for  the  stage  was 
approximately  4.12  and  60/40  respectively.  The  high  rotor 
solidity  was  necessary  to  limit  the  DF  requirement  of  the  rotor 
and  even  with  this  high  a  solidity,  the  rotor  DF  was 
approximately  0.7. 

As  shown  earlier.  Figure  8  gives  the  non-dimensional 
power  loss  of  the  secondary  stream  for  the  flow  control  process. 
An  adiabatic  efficiency  may  be  assigned  to  this  value  to 
account  for  losses  within  a  “hypothetical'’  secondary  flow 
delivery  system  as  shown  in  Figure  2  and  give  an  “actual” 
power  input  for  the  entire  secondary  flow  system,  yielding: 


j,  actual 


Vi 


(9) 


Combining  this  with  the  power  input  of  the  rotor  and  using 
the  definition  of  adiabatic  efficiency  gives  the  total  power 
requirement  of  the  stage  shown  in  Figure  2  as: 

(  r~ 1 

r  -1 


71  „ 


P  =  " \CpTt,c,t 


V 


J 


+  tKCPTt,c,  2' 


(10) 


Vc  Vj 

The  first  term  in  Eq.  (10)  is  the  power  input  of  the  core 


flow  from  inlet  (c,l)  to  exit  (c,3).  The  core  pressure  ratio,  Jt c , 

includes  the  “apparent”  loss  through  the  flow  controlled  stator. 
Since  this  flow  control  technique  imparts  momentum  to  the 
core  flow,  the  “apparent”  total  pressure  loss  through  the  stator 
is  negligible.  The  CFD  simulations  of  the  stator  section 
showed  a  reduction  of  only  0.27%  in  the  core  mass  averaged 
total  pressure  for  the  momentum  coefficient  of  0.1456. 


Therefore,  the  core  total  pressure  ratio,  71  ,  could  be 
considered  equivalent  to  the  rotor  pressure  ratio. 


The  core  efficiency,  77  ,  is  another  matter.  This  efficiency 

could  be  affected  by  heat  transfer  from  the  secondary  flow 
stream  to  the  core  flow  making  the  process  non-adiabatic 
through  the  stator.  This  violates  the  assumption  made  in 


deriving  the  overall  power  input  equation  which  used  the 
definition  of  adiabatic  efficiency  for  the  core  flow  power  input. 
Also,  the  less  efficient  the  secondary  flow  system,  i.e.  the 

smaller  Tj  . ,  the  larger  the  secondary  flow  injection  temperature 

for  a  given  secondary  flow  total  pressure  ratio.  This  would 
make  the  core  flow  efficiency  in  this  analysis  a  function  of  the 
secondary  flow  efficiency  due  to  heat  transfer  effects.  An 
assessment  of  this  effect  at  this  time  is  not  accounted  for  in  this 
analysis  since  the  CFD  simulations  presented  earlier  were 
performed  with  the  injection  total  temperature  equal  to  the  core 
flow  total  temperature.  For  this  analysis,  this  efficiency  will  be 
held  at  the  rotor  efficiency  value  of  the  meanline  analysis,  i.e. 


VC=V 


rotor  * 


Continuing  with  the  analysis,  the  adiabatic,  isentropic 
power  input  of  the  stage  is  the  first  term  of  Eq.  (10)  with  the 
adiabatic  efficiency  set  to  one: 

(  y- 1  b 

7  -1 


P  =  iji  r  T 

1  isen 


7tn 


V 


(11) 


J 


The  ratio  of  Eq.  (11),  the  isentropic  power,  to  Eq.  (10),  the 
actual  power  input,  defines  an  efficiency  and  is,  upon 
substituting  and  re-arranging: 

P  77 

V  =  ^T  =  7 - 7 - - - 7 - A  (12) 


\  +  P. 


1 

r- 1 


V^C 
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-1 


1 

Vc 


\  7 

f  \ 


V  V  j  J 
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Figure  10  shows  Eq.  (12)  plotted  against  the  secondary 

_  * 

flow  system  efficiency,  77  . .  The  value  for  Pj  jsen  was  chosen 


from  Figure  8  for  a  DF  of  0.98.  Two  values  of  overall 
efficiency  are  highlighted  in  Figure  10.  There  is  no  particular 
significance  for  these  two  values,  they  were  chosen  for 
comparison  purposes  only  to  assess  the  sensitivity  of  the  overall 
efficiency  to  the  secondary  flow  efficiency  raised  in  private 
discussions  [17].  As  can  be  seen,  the  overall  efficiency  is  0.883 
for  a  secondary  flow  system  efficiency  of  0.7.  There  is  1.7 
points  in  overall  efficiency  taken  up  in  powering  the  secondary 
flow  system.  A  28.5%  reduction  in  the  secondary  flow 

efficiency,  i.e.  7 J  ■  =  0.50 ,  results  in  an  overall  efficiency  of 


0.8764,  a  reduction  of  0.75%  over  the  Tj  -  =0.70  case.  The 

sensitivity  of  the  overall  efficiency  to  the  secondaiy  flow 
system  efficiency  is  not  very  large  in  the 

0.5  <  IJj  <0.70  range  for  the  secondary  flow  power  input 

under  consideration  and  for  this  stage  configuration.  From  this 
simple  analysis,  the  penalty  to  the  stage  using  this  flow  control 
technique  does  not  appear  to  have  a  severe  impact  on 

performance,  even  with  low  secondary  flow  system  efficiency. 
Although,  as  mentioned  earlier,  this  does  not  take  into  account 
the  non-adiabatic  effects  due  to  heat  transfer  between  the 
secondary  flow  and  core  flow. 
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Figure  10  Overall  Efficiency  vs  Secondary  Flow  System 
Efficiency 


SYSTEM  LEVEL  ANALYSIS 

The  overall  engine  performance  implication  of  utilizing  the 
flow  control  technique  for  front  fan  applications  was 
investigated  as  part  of  this  effort.  A  two-spool,  mixed-flow, 
augmented  turbofan  with  a  year-2000  production-capable  level 
of  technology  is  used  for  analysis  purposes.  Table  3  provides 
the  basic  thermodynamic  cycle  characteristics  of  this  engine, 
and  Figure  1 1  provides  a  graphical  representation  of  the 
geometry/flowpath.  (It  is  important  to  note  that  this 
configuration  is  notional,  and  does  not  represent  a  fielded 
engine  or  an  engine  in  development.) 


Table  3  Notional  Engine  Cycle  Characteristics 


OPR 

20 

BPR 

0.27 

m 

124.7  kg/sec 

Throttle  Ratio 

1.15 

Inlet  Recovery 

0.97 

n  fan 

4.44 

V  fan 

0.87 

Vhpc 

4.51 

Whpc 

0.86 

Maximum  T3 

866.5  K 

Maximum  T4i 

1894.3K 

n hpt 

2.29 

Vhpt 

0.88 

f hpt , vane ,  cooling 

10.00%  of  core  flow 

f hpt ,  blade ,  cooling 

5.50%  of  core  flow 

n iPt 

2.03 

V,P< 

0.89 

f Ipt, vane,  cooling 

1 .60%  of  core  flow 

f Ipt, blade, cooling 

1 .40%  of  core  flow 

Nozzle  CFG 

0.95 

FnSLS-mil 

98875  N 

FfrSLS-max 

147841 N 

SFCsLS-mil 

0.899  kg/hr/daN 

Figure  11  Notional  Engine  Cross  Section 


Table  4  describes  the  pertinent  aerodynamic  data  related  to 
the  3-stage  front  fan  design  in  the  notional  engine.  Of 
particular  relevance  to  the  current  study  is  the  average  amount 
of  stator  turning  at  the  hub  locations. 


Table  4  Fan  Aerodynamic  Parameters  of  Notional  Engine 


Corrected  Flow  /  Annulus  Area 

170.9  kg/sec/mA2 

Number  of  Stages 

3 

Stage  Loading  (2gJ(delH)/UA2) 

0.83 

Inlet  Corrected  Tip  Speed 

453.8  m/s 

Inlet  Physical  Tip  Speed 

454.8  m/s 

Exit  Hub  Speed 

324.9  m/s 

Average  Hub  Stator  Turning 

56.6  deg 

In  general,  without  the  use  of  active  flow  control, 
separation  losses  limit  the  amount  of  stator  turning  at  the  hub 
location  of  any  stage  to  -65  degrees.  Through  the  use  of  active 
flow  control,  this  constraint  becomes  less  severe  and  fan 
aerodynamic  designs  with  reduced  stage  count  become  feasible. 
The  weight  savings  associated  with  a  stage  reduction  are 
countered  by  thermodynamic  losses  associated  with  powering 
the  flow  control  system.  The  magnitudes  of  these 
savings/losses  are  quantified  as  follows. 

The  compressor  map  of  the  fan  itself,  i.e.  all  stages,  is 
assumed  to  be  unchanged,  and  horsepower  is  extracted  from  the 
low-pressure  turbine  to  power  the  flow  control  system.  Figure 
12  shows  the  effect  of  low-spool  horsepower  extraction  on 
thrust  and  specific  fuel  consumption  (each  at  sea-level  static 
conditions  without  the  use  of  an  afterburner). 
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Figure  12  Thrust  and  SFC  versus  Secondary  Flow  Power 
Extraction 


A  reasonable  assumption  for  the  power  required  to  drive 
the  flow  control  system  is  5%  of  the  fan  power  in  the  notional 
engine.  This  5%  requirement  was  arrived  at  by  examining  the 
ratio  of  the  secondary  flow  power  input  (second  term  of  Eq. 
(10))  to  core  fan  power  requirement  (first  term  of  Eq.  (10)) 
from  the  ID  power  analysis  outlined  in  the  previous  section: 
f  \ 
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Using  the  fan  efficiency  and  pressure  ratio  from  Table  3 
along  with  setting  the  secondary  flow  efficiency  to  0.5  and 

Pj  jsen  to  0.0048  yields  the  approximate  value  of  0.025  for 


Multiplying  this  by  two  for  additional  losses  (e.g. 


bearing,  shaft,  etc.)  gives  a  conservative  estimate  of  0.05  or  5%. 

Will  all  this  in  mind,  this  yields  a  requirement  of 
approximately  1 120  kW,  and  results  in  a  (dry)  net  thrust  loss  of 
1463  N  and  specific  fuel  consumption  (SFC)  increase  of  0.013 
kg/hr/daN.  (Note  that  at  sea  level  static  in  full  afterburner  with 
1120  kW  extraction,  the  net  thrust  loss  is  1730  N  when 
compared  to  the  notional  engine  in  full  afterburner.) 

A  reduction  in  the  number  of  fan  stages  from  3  to  2  will 
result  in  a  fan  weight  savings  in  the  1112-1245  N  regime. 
Figure  13  shows  the  influence  of  average  stage  loading 
coefficient  on  1st  stage  stator  turning  angle  at  the  hub  location. 
For  loadings  above  0.95,  Figure  13  shows  that  a  two  stage 
design  will  require  an  amount  of  hub  stator  turning  greater  than 
65  degrees.  The  relative  Mach  number  at  the  first  stage  rotor 
tip  is  shown  for  the  2  and  3  stage  designs  as  a  function  of 
loading  coefficient  in  Figure  14.  The  average  stage  loading 
coefficient  must  increase  from  0.8  to  1.4  in  order  for  the  2  stage 
design  to  have  the  same  1st  stage  rotor  tip  relative  Mach  number 
as  the  3  stage  design.  An  increase  in  average  stage  loading 
coefficient  implies  a  higher  stator  diffusion  requirement,  an 
environment  in  which  flow  control  can  help. 

The  system  level  analysis  suggests  that  active  fan  flow 
control  can  be  used  to  increase  overall  engine  thrust  to  weight. 
The  notional  engine  had  a  thrust-to-weight  ratio  of  5.62.  A 


replacement  of  the  3  stage  front  fan  with  a  2  stage  design  that 
incorporates  active  flow  control  will,  despite  the 
thermodynamic  losses,  yield  an  overall  engine  thrust-to-weight 
improvement  to  the  5.80  regime.  Further  analysis  is  required  to 
quantify  the  benefit  for  other  engine  cycles,  particularly  ones 
that  involve  fixing  the  engine  volumetric  size  and  allowing  for 
higher  OPR.  It  is  in  these  cycles  that  one  may  see  a  SFC 
benefit  while  maintaining  or  slightly  increasing  the  thrust  to 
weight  ratio. 


Loading 


Figure  13  First  Stage  Stator  Hub  Turning  versus  Loading 


Figure  14  Relative  Mach  Number  versus  Loading 

EXPERIMENTAL  RESULTS 

Cascade 

A  cascade  was  created  by  simply  extruding  the  cross 
section  outlined  earlier  with  internal  cavities  sufficiently  sized 
to  provide  the  necessary  mass  flow  to  feed  the  injection  and 
suction  slots  on  the  blade  surface.  The  blade  section  with  the 
larger  injection  slot  height  to  pitch  ratio  (h/s)  of  .011974  was 
chosen.  The  cross  section  of  the  cascade  is  shown  in  Figure  15 
and  the  solid  model  in  Figure  16.  Secondary  flow  to  the 
cascade  is  provided  from  only  one  side  of  the  blade  to  allow 
viewing  from  the  other  side  of  the  cascade  test  section  for 
DPIV  measurements.  Also,  this  would  be  the  intended 
configuration  for  a  stage  as  shown  in  Figure  2.  Although,  the 
blade  cross  section  provided  sufficient  area  to  accommodate  the 
necessary  mass  flow,  it  was  equal  to  the  respective  slot  areas  to 
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which  they  supplied  air  flow.  Ideally  one  would  want  the  cross 
sectional  areas  supplying  the  secondary  flow  to  be  larger,  but 
that  was  not  possible  due  to  geometric  constraints  and  structural 
rigidity  considerations.  Although  manufacturing  the  cascade 
using  EDM  technology  is  a  fairly  easy  task,  a  minimum  wall 
thickness  is  required  to  prevent  warping  and  maintain 
tolerances. 

The  blowing  and  suction  supply  cavities  need  to  be  located 
within  the  central  portion  of  the  blade  cross  section  since  this 
portion  of  the  blade  provides  sufficient  metal  thickness.  With 
the  blowing  and  suction  slots  located  to  the  ends  of  the  blade, 
this  requires  the  internal  flow  to  negotiate  an  undesirable  180 
degree  turn  as  can  be  seen  in  Figure  15.  Unfortunately,  without 
very  careful  design  considerations  in  this  region,  this  will  result 
in  significantly  more  loss  generation.  Insufficient  time  did  not 
allow  for  any  detailed  analysis  of  this  area.  With  all  the 
geometric,  manufacturing  and  time  constraints,  the  cross 
section  in  Figure  15  is  considered  the  best  that  could  be  done. 

The  cascade  test  section  is  six  inches  tall  and  contains  six 
blades  across  the  test  section  with  a  blade  pitch  of  0.0762 
meters.  The  blade  chord  measured  as  a  straight  line  from 
leading  to  trailing  edge  is  0.127  meters  giving  a  solidity  of 
1.67.  The  cascade  inlet  Mach  number  was  0.7  with  an  inlet 
angle  of  68  degrees. 


Suction  Slot 


Figure  16  Flow  Control  Cascade 


Secondary  Flow  System 

General  Layout 

The  secondary  flow  system  is  created  using  a  centrifugal 
automotive  supercharger  powered  by  an  electric  motor  and 
controlled  by  an  AC  inverter.  Figure  17  gives  a  schematic  of 
this  system  along  with  measurements  made  at  various  stations 
along  the  flowpath  of  the  system.  The  by-pass  valve  shown 
will  be  explained  later. 


Figure  17  Experimental  Secondary  Flow  Circuit 

The  supercharger  was  chosen  because  it  could  provide  the 
desired  pressure  ratio  and  power  input  for  the  momentum 
coefficients  of  interest.  It  was  difficult  to  locate  a  supercharger 
with  sufficient  surge  margin  for  the  flow  rate  of  interest.  The 
cascade  corrected  flow  rate  is  approximately  4.80  kg/sec.  Since 
only  4  blades  out  of  the  6  (i.e.  5  passages)  are  controlled,  the 
flow  rate  that  needs  to  be  controlled  is  3.84  kg/sec.  The  desired 
fraction  of  secondary  flow  to  core  flow  (i.e.  flow  fraction,  f) 
from  Table  1  is  0.0454,  which  means  that  the  required 
corrected  flow  rate  for  the  secondary  flow  system  is  0.174 
kg/sec.  Searching  all  the  available  centrifugal  supercharger 
maps  puts  this  flow  rate  very  close  to  the  surge  line  for  the  even 
the  smallest  supercharger  that  could  provide  the  necessary 
pressure  ratio.  A  detailed  map  of  the  supercharger  chosen  for 
this  experiment  was  not  available,  but  a  map  of  a  “very  close” 
competitor  with  “similar  performance”  was  obtained  and 
verified  that  it  would  be  operating  very  close  to  surge  at  this 
flow  rate. 

The  desired  pressure  ratio  of  the  secondary  flow  system 
from  Table  1  is  approximately  1.48.  This  pressure  ratio  is  the 
ratio  of  the  blowing  jet  flow  total  pressure  to  the  core  flow  total 
pressure.  In  actuality,  the  secondary  flow  system  pressure  ratio 
is  equal  to  the  ratio  of  the  suction  slot  total  pressure  to  the  inlet 
blowing  slot  total  pressure.  This  pressure  ratio  is 
approximately  1.4.  But  this  pressure  ratio  does  not  take  into 
account  any  flow  losses  within  the  secondary  flow  system. 

Due  to  short  period  of  time  afforded  for  the  design  of  the 
secondary  flow  system,  a  simple  dump  manifold  was  created  to 
supply  the  air  to  the  cascade  section.  Figure  18  shows  a  picture 
of  the  dump  manifold  with  the  cover  removed.  Of  course  this 
type  of  manifold  is  not  desirable  from  a  loss  perspective,  so  an 
analysis  was  conducted  to  determine  the  required  pressure  ratio 
of  the  centrifugal  compressor  needed  to  deliver  a  1.4  pressure 
ratio  on  the  cascade  blade.  The  loss  analysis  showed  that  the 
centrifugal  compressor  needed  to  supply  a  pressure  ratio  on  the 
order  1.9  to  achieve  the  desired  blowing  total  pressure.  This 
places  the  compressor  in  surge  for  a  corrected  flow  rate  of 
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0.174  kg/sec.  Therefore,  a  by-pass  valve  was  added  to  the 
system  as  shown  in  Figure  17  to  re-circulate  flow  and  keep  the 
centrifugal  compressor  from  surging  while  providing  sufficient 
pressure  ratio  to  the  secondary  flow.  This  assumes  that  the 
inlet  blowing  slots  for  the  flow  control  are  choked  and  that  re¬ 
heating  effects  of  the  re-circulating  secondary  flow  for  the  short 
blow  down  test  time  of  approximately  10  seconds  would  not 
have  a  significant  effect. 

Drive  System 

The  approximate  amount  of  power  required  for  the 
secondary  flow  system  can  be  calculated  from  a  simple  ID 
analysis  as  carried  out  in  the  previous  sections  and  re-written 
here: 

f  r- 1  A 

n  7  -1 

P  =  mcTt -  (14) 

7 

The  blow-down  tunnel  is  a  pressurized  system.  The  total 
pressure  is  anywhere  from  1. 378-1. 723x10s  Pa  from  past 
experience  with  the  specified  inlet  condition  to  this  cascade. 
The  tunnel  total  temperature  can  be  approximated  as  a  standard 
day  condition.  With  this  tunnel  total  pressure  and  temperature 
and  a  corrected  mass  flow  of  0.174  kg/sec,  the  actual  mass  flow 
in  the  secondary  system  may  vary  between  0.238  kg/sec  to 
0.295  kg/sec.  Approximating  the  compression  efficiency  to  be 
70%  from  a  comparable  compressor  map  would  lead  to  a  power 
requirement  of  18.64  to  23.86  kW.  A  3400  RPM,  37.28  kW, 
460  Volt,  3  phase  electric  motor  was  chosen  to  be  sufficient  for 
the  experiment  controlled  by  an  AC  inverter  capable  of 
continuous  44.74  kW  control.  . 

The  power  transmission  to  the  supercharger  was 
accomplished  by  using  a  high  precision  timing  pulley  and  belt 
arrangement.  The  supercharger  speed  needed  to  be 
approximately  42,000  rpm  in  order  to  achieve  the  required 
pressure  ratio.  There  is  an  internal  step  up  within  the 
supercharger  of  4.1:1  thereby  requiring  an  external  pulley 
arrangement  of  3:1  using  a  maximum  motor  speed  of  3400 
RPM.  At  these  speeds,  stock  steel  pulleys  are  close  to  their 
maximum  rim  speed,  but  a  system  was  found  that  could 
transmit  the  required  power  and  was  used  successfully  in  the 
experiment. 

DPIV  Results 

The  desired  inlet  conditions  for  the  cascade  test  were  an 
inlet  angle  of  68  degrees  and  Mach  number  of  0.7. 
Unfortunately,  the  secondary  flow  system  choked  at  less  than 
half  the  desired  flow  rate  of  approximately  0.238  kg/sec.  A 
simple  test  was  performed  to  determine  which  part  of  the 
secondary  flow  system  was  limiting  the  flow  rate.  Two  of  the 
blade  surface  injection  slots  were  taped  off  and  the  flow  rate 
remained  unchanged.  The  tape  was  removed  and  two  of  the 
blade  surface  suction  slots  were  taped  off  which  in  turn  halved 
the  flow  rate.  Therefore,  from  this  simple  analysis  it  was 
determined  that  the  somewhere  within  the  suction  side  of  the 
secondary  flow  system  choking  occurred.  Currently  a  CFD 
analysis  of  the  blade  passages  is  being  performed  to  determine 
where  within  the  suction  portion  of  the  secondary  flow  system 
this  occurred. 


Figure  18  Dump  Manifold  Picture 


With  the  limited  amount  of  flow  rate  in  the  secondary 
system,  the  inlet  Mach  number  to  the  cascade  was  lowered  to 
0.3  in  order  to  try  and  match  the  flow  fraction  and  momentum 
coefficient  of  0.454  and  0.197  respectively  of  the  CFD 
simulations  outlined  earlier.  This  was  the  lowest  Mach  number 
the  cascade  tunnel  could  run  stably.  Also,  to  achieve  the 
desired  pressure  ratio  of  1 .87  within  the  secondary  flow  system, 
the  by-pass  valve  was  slightly  opened  to  re-circulate  a  fraction 
of  the  flow  prior  to  the  manifold  to  keep  the  supercharger  from 
surging.  An  attempt  was  made  to  calculate  a  momentum 
coefficient  by  measuring  the  static  pressure  at  the  endwalls  on 
each  side  of  the  blade  injection  slots.  One  could  then  calculate 
a  momentum  coefficient  assuming  the  total  pressure  and 
temperature  measurement  within  the  blowing  plenum  portion  of 
the  dump  manifold  shown  in  Figure  17  (location  3  in  the  figure) 
to  be  the  total  pressure  and  temperature  at  the  injection  slot. 
Unfortunately,  the  thermocouple  measurements  used  to 
determine  the  total  temperature  at  the  inlet  to  and  within  the 
blowing  plenum  of  the  dump  manifold  failed  during  the  test. 

Data  taken  for  the  secondary  flow  system  prior  to  the  test 
showed  that  the  supercharger  could  operate  in  the  by-pass 
configuration  at  an  approximate  adiabatic  efficiency  of  0.65  for 
a  total  pressure  ratio  of  1.8.  This  was  the  operating  condition 
chosen  for  the  test.  The  measured  inlet  temperature  to  the 
supercharger  at  this  operating  condition  during  the  test  was  309 
K.  Assuming  an  adiabatic  efficiency  of  0.65  with  a  total 
pressure  ratio  of  1.8  puts  the  exit  temperature  of  the 
supercharger  at  396  K.  With  this  temperature,  along  with  the 
measured  total  and  static  pressure  prior  to  the  blowing  plenum 
(location  2  in  Figure  17)  to  calculate  the  mass  flow  of  the 
secondary  system,  yielded  a  secondary  mass  flow  rate  of 
0.0716  kg/sec.  The  uncertainty  in  this  calculation  was 
determined  to  be  +0.001738  kg/sec.  The  core  flow  rate  for  the 
blow  down  tunnel  operating  at  a  Mach  number  0.3  was 
calculated  to  be  2.0  kg/sec.  This  yields  a  flow  fraction  of 
secondary  to  primary  flow  of  0.0358.  This  is  approximately 
1%  of  the  core  flow  less  than  the  desired  value  of  0.0454.  It 
was  also  noted  that  measurements  taken  prior  to  the  failure  of 
the  thermocouples  showed  a  significant  heat  loss  in  the  dump 
manifold.  That  combined  with  the  failure  of  the  temperature 
measurements  during  the  test  made  any  attempt  to  calculate  the 
momentum  coefficient  unreliable. 

Figure  19  shows  the  interrogation  regions  for  the 
DPIV  measurements.  The  interested  reader  is  pointed  to 
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references  [2,18]  for  a  more  detailed  look  at  the  DPIV  setup. 
Only  regions  A,  B  and  C  were  obtained  due  to  time  constraints. 

There  are  many  factors  associated  with  the  DPIV 
uncertainty-calculation  process  (laser,  CCD,  seeding,  imaging, 
algorithms,  oscilloscope,  etc).  The  highest  uncertainty  was 
found  to  be  associated  with  the  velocity  calculation  which 
involves  Ax  (the  displacement  in  pixels  of  each  interrogation 
region).  At  (the  time  interval  between  the  two  exposures),  and 
the  magnification  of  the  digital  image  relative  to  the  object 
(pix/m).  The  displacement  in  pixels  obtained  by  peak-locator 
algorithms  can  provide  sub-pixel  accuracy  (<  0.1  pixels)  after 
correction  for  various  biases.  The  At  was  adjusted  to  yield 
typical  displacements  of  the  main  stream  >10  pixels,  and  the 
uncertainty  is  thus  <1%.  Values  in  the  wake  region,  however, 
may  have  higher  uncertainties  due  to  the  lower  Ax.  The 
maximum  uncertainty  in  the  At  was  calculated  from  the  time 
interval  between  the  two  laser  pulses  with  the  aid  of  an 
oscilloscope  (uncertainty  2%).  It  was  found  that  this 
uncertainty  increases  with  lower  laser  power  and  with  lower  At. 
A  conservative  number  for  the  present  experiments,  which 
employed  a  At  of  about  4  us  and  powers  around  10  mJ,  was 
found  to  be  <1%.  The  magnification  was  measured  using 
images  of  grids  located  in  the  laser-sheet  plane  to  better  than 
1%.  Combining  these  conservative  measurements  of 
uncertainty  yields  a  maximum  error  of  <  2%  for  the  free-stream 
regions  and  -10%  in  the  low  speed  areas  such  as  the  wake 
region. 


A  composite  image  of  the  interrogated  regions  for  the 
baseline  case  (i.e.  no  flow  control)  is  shown  in  Figure  19. 


As  you  can  see,  a  significant  amount  of  separation  occurs 
along  the  suction  surface  of  the  blade.  Similarly,  a  composite 
image  of  the  controlled  blade  is  shown  in  Figure  20.  Again, 
separation  is  seen,  but  at  a  significantly  lower  level  over  the 
baseline  case,  showing  that  the  control  is  having  a  desirable 
effect. 

Significant  difference  were  seen  in  the  injection  slot  static 
pressure  measurements  taken  at  each  endwall  which  leads  to 
the  conclusion  of  very  non-uniform  injection  and  a  high 
probability  for  non-uniform  suction.  A  detailed  3D  CFD 
analysis  is  underway  to  model  this  cascade  section  and 
preliminary  results  are  confirming  this  hypothesis.  Any  radial 


injection  would  imply  a  reduction  in  momentum  coefficient 
aligned  with  the  core  flow.  A  subsequent  redesign  which  will 
seek  to  provide  injection  and  suction  uniformity  will  follow  the 
analysis. 


Figure  20  DPIV  Experimental  Baseline  Composite  DPIV 
Image 


Figure  21  Experimental  Flow  Control  Composite  DPIV 
Image 


CONCLUSIONS 

A  continuous  steady  flow  control  concept  using  an  inset 
cavity  and  pumping  system  to  continually  re-circulate  a 
secondary  flow  which  in  turn  enhances  the  core  flow  diffusion 
was  presented.  A  control  volume  analysis  was  performed  to 
find  relevant  non-dimensional  parameters  and  a  CFD 
simulation  executed  varying  those  parameters  showed  that  an 
apparent  DF  of  0.98  could  be  achieved.  The  word  “apparent”  is 
used  because  the  secondary  flow  control  stream  adds 
momentum  to  the  core  flow  with  the  secondary  flow  system 
requiring  work  input  through  a  pumping  system.  The  “penalty” 
for  the  secondary  system  is  addressed  with  a  simple  power  and 
availability  analysis.  The  availability  analysis  showed  that  the 
process  involved  with  this  technique  resulted  in  approximately 
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a  50%  lower  power  loss  through  the  stator  over  the  separated 
baseline  case.  This  analysis  included  the  losses  within  the 
secondary  flow  stream,  but  did  not  account  for  losses  within  a 
delivery  system  of  this  secondary  flow.  This  means  that 
thermodynamically  this  flow  control  process  is  more  efficient 
than  the  separated  baseline  case,  but  losses  within  a  delivery 
system  could  make  the  two  thermodynamically  equivalent. 
Even  with  thermodynamic  equivalence,  the  flow  controlled 
case  still  has  the  advantage  of  achieving  the  desired  diffusion 
levels  and  exit  conditions  over  the  uncontrolled  baseline.  The 
ID  power  analysis  was  used  to  assess  this  technique’s  impact 
on  the  efficiency  of  an  axial  compression  stage  and  the 
sensitivity  of  this  efficiency  to  the  secondary  flow  system’s 
efficiency.  A  stage  with  a  pressure  ratio  of  4.12  and  a 
compressor  adiabatic  efficiency  of  0.90,  yielded  an  overall 
efficiency  of  0.883  for  a  secondary  flow  system  efficiency  of 
0.7.  Lowering  the  secondary  flow  system  efficiency  to  0.5 
gave  a  0.876  overall  efficiency.  The  sensitivity  of  the  overall 
efficiency  to  the  secondary  flow  system  efficiency  was  not 
excessively  large.  Therefore,  it  may  be  feasible  to  achieve  a 
4.0  pressure  ratio  with  a  0.88  adiabatic  efficiency. 

Also,  a  system  level  analysis  is  presented  to  assess  the 
merits  that  may  be  realized  in  a  notional  engine  with  this  type 
of  flow  control.  The  notional  engine  contained  a  3  stage  fan 
which  was  replaced  by  a  2  stage  fan  using  flow  control.  The 
power  to  drive  the  flow  control  system  was  taken  off  the  low 
pressure  turbine.  For  this  preliminary  study,  an  engine  fuel 
consumption  penalty  and  a  thmst-to-weight  benefit  were 
realized  with  a  specific  fuel  consumption  and  thrust-to-weight 
ratio  increase  of  1.5%  and  3.2%  respectively. 

A  cascade  experiment  was  performed  to  demonstrate  the 
concept  and  was  conducted  in  a  blow-down  cascade  tunnel. 
The  secondary  flow  system  was  designed  using  an  automotive 
supercharger  driven  by  an  electric  motor  to  provide  continuous 
re-circulating  flow.  The  experiment  was  not  successful  in 
providing  the  simulated  levels  of  diffusion  due  to  flow 
limitations,  i.e.  choking,  within  the  suction  portion  of  the  flow 
control  scheme.  Even  with  this  limitation,  significant 
improvements  in  diffusion  were  qualitatively  seen  from  the 
DPIV  measurements.  It  is  believed  that  non-uniform  injection 
and  suction  led  to  the  less  than  expected  behavior  of  the  flow 
control  technique.  A  3D  CFD  analysis  is  ongoing  and 
preliminary  findings  are  showing  large  injection  and  suction 
non-uniformities.  A  redesigned  cascade  section  with  improved 
injection  and  suction  characteristics  will  follow. 
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NOMENCLATURE 


AR 

aspect  ratio 

AVR 

axial  velocity  ratio 

BPR 

bypass  ratio 

c. 

injection  momentum  coefficient 

Cr 

static  pressure  rise  coefficient 

DF 

diffusion  factor 

/ 

flow  fraction 

h/S 

ratio  of  injection  slot  height  to  pitch 

F 

integrated  loading 

Fn 

thrust 

H, 

stagnation  enthalpy 

m 

mass  flow 

m 

A 

specific  mass  flow 

V 

velocity  magnitude 

OPR 

overall  pressure  ratio 

P 

power 

S 

pitch 

s 

entropy 

SFC 

specific  fuel  consumption 

aM 

mass  averaged  flow  angle 

r 

circulation 

?7 

adiabatic  efficiency 

n 

total  pressure  ratio 

cr 

solidity 

P 

density 

T 

total  temperature  ratio 

4* 

availability 

Subscripts 


c 

core  flow 

.i 

secondary  flow 

X 

axial  direction 

y 

tangential  direction 

t 

stagnation  condition 

fan 

Fan 

hpc 

high  pressure  compressor 

hpt 

high  pressure  turbine 

lpt 

low  pressure  turbine 

sis 

sea  level  static  conditions 
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